Simulations and Experiments on Polarisation Squeezing in Optical Fibre 
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We investigate polarisation squeezing of ultrashort pulses in optical fibre, over a wide range of 
input energies and fibre lengths. Comparisons are made between experimental data and quantum 
dynamical simulations, to find good quantitative agreement. The numerical calculations, performed 
using both truncated Wigner and exact -\-P phase-space methods, include nonlinear and stochastic 
Raman effects, through coupling to phonons variables. The simulations reveal that excess phase 
noise, such as from depolarising GAWBS, affects squeezing at low input energies, while Raman 
effects cause a marked deterioration of squeezing at higher energies and longer fibre lengths. The 
optimum fibre length for maximum squeezing is also calculated. 
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I. INTRODUCTION 

The search for efficient means of quantum squeezing, 
in which quantum fluctuations in one observable are re- 
duced below the standard quantum limit, at the expense 
of increased fluctuations in the conjugate, has been at 
the heart of modern developments in quantum optics [||. 
As well as for the fundamental interest of highly nonclas- 
sical light, optical squeezing is of interest for quantum 
information applications. Possible uses include: generat- 
ing entanglement for quantum communication [3], making 
measurements below the standard quantum limit [3|, and 
for precise engineering of the quantum states of matter [J| . 

The use of optical fibre for quantum squeezing has 
considerable technological advantages, such as generat- 
ing squeezing directly at the communications wavelength 
and use of existing transmission technology. There is, 
however, a significant disadvantage in the excess phase 
noise that arises from acoustic waves, molecular vibra- 
tions, and defects in the amorphous silica. 

Here we present an in-depth numerical and experimen- 
tal study of polarisation squeezing in a single-pass scheme 
that successfully reduces the impact of this excess phase 
noise. The numerical simulations represent a quantita- 
tive, experimentally testable solution of quantum many- 
body dynamics. 

The first proposals for the generation of squeezed light 
using the x^^"* nonlinearity date back to 1979, with 
schemes involving a nonlinear Kerr interferometer [5[ 
or degenerate four- wave mixing [6|. The first experi- 
mental demonstration used four-wave mixing in atomic 
samples Q- The Kerr effect in optical fibers was also 
proposed as a mechanism for squeezing light @, M> [l3l ■ 
Squeezing using fibres was first successfully implemented 
using a continuous wave laser, and was observed by a 



phase shifting cavity |ll| . 

However, early experiments |8l. IqI. I 111] were severely lim- 
ited by the phase noise intrinsic to optical fibre. Such 
noise occurs in the form of thermally excited refractive 
index fluctuations in the fiber 'l_^, ll3[, and arises from 
Guided Acoustic Wave Brillouin Scattering (GAWBS) 
and 1// noise. A substantial theoretical breakthrough 
was the recognition that short pulses - ideally in the form 
of solitons - could lead to much higher peak powers, thus 
allowing the generation of nonclassical light in with fiber 
lengths short enough so that thermally induced phase 
noise was not an issue. Such short pulses required a 
true multi-mode theoretical approach[14j, which led to 
the first predictions of pulsed squeezing, and to an un- 
derstanding of the scaling laws involved [l^ . 

These predictions were confirmed in a landmark ex- 
periment by Rosenbluh and Shelby [l6| , which used in- 
tense, sub-picosecond laser pulses to eliminate much of 
the phase noise, and a simpler interferometric setup [lO| 
in a balanced configuration. All fiber squeezers since have 
exploited ultrashort pulses. Observation schemes imple- 



mented with standard fibers include 
cavities [ll|, ii) spectral filtering [l 
iii) balanced interferometer s Im |2 
asymmetric interferometers 1261 1271 . |2 
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two-pulse, single-pass method generatim 
uum [3l|, [34I or polarization squeezing [3: 
Squeezing the polarization variables 
promising alternative 35| to the squeezing in the ampli- 
tude quadrature or the photon number, which the vast 
majority of fiber squeezing experiments until now have 
implemented. That the quantum polarization variables 
could also exhibit noise reduction was first suggested by 
Chirkin et al. in 1993 [36[. This proposal for polar- 
ization squeezing is similar to earlier remarks by Walls 
and ZoUer concerning atomic spin squeezing '37| due to 
the systems' mathematical similarities. The first exper- 
iment to exploit the quantum properties of polarization 
was performed by Grangier et al. in 1987 in a squeezed- 



light-enhanced polarization interferometer [3^| . The first 
explicit demonstration was achieve d by S0rensen et al. 
in the context of quantum memory [39| . Such a promis- 
ing application sparked intensified interest, resulting in 
a number of theoretical investigations, e.g. ^35 . i40l . |4l|. 
These in turn spawned a plethora of experiments in a 
variety of different systems: optical parametric oscilla- 
tors [42, 43, 44], optical fibers [29, 33, 34] and cold atomic 
samples [45]. 

In this paper we present a detailed experimental and 
theoretical investigation of the single-pass method for 
creating polarization squeezing, building upon our pre- 
vious work |33l . 1461 ]. This efficient and novel squeezing 
source has a number of advantages compared with previ- 
ous experiments producing bright squeezing. For exam- 
ple, this setup is capable of producing squeezing over a 
wide range of powers, in contrast to asymmetric Sagnac 
loop schemes. There is thus a certain similarity to exper- 
iments using a Mach Zehnder interferometer as a flexi- 
ble asymmetric Sagnac loop [28|. The interference of a 
strong squeezed and a weak 'coherent' beam in asymmet- 
ric loops however gives rise to a degradation in squeezing 
due to the dissimilarity of the pulses as well as losses from 
the asymmetric beam splitter. 

In the single-pass scheme, the destructive effect aris- 
ing from interfering dissimilar pulses (in power, temporal 
and spectral shape) is avoided by interfering two strong 
Kerr-squeezed pulses that co-propagate on orthogonal 
polarization axes. For equal power they have been found 
to be virtually identical within measurement uncertain- 
ties in, e.g. spectrum, autocorrelation and squeezing. 
This scheme presents the potential to measure greater 
squeezing and provides a greater robustness against input 
power fluctuations. Formally this interference of equally 
squeezed pulses is reminiscent of earlier experiments pro- 
ducing vacuum squeezing, for example 16, 22, 32]. The 
advantage here is that no extra local oscillator is needed 
in the measurement of polarization squeezing. 

These novel experiments allow a careful experimental 
test of the multi-mode theory of optical squeezing. Here 
we make use of the comprehensive model developed by 
Carter and Drummond[47J that includes the electronic 
X^^' nonlinear responses of the material and nonresonant 
coupling to phonons in the silica. The phonons provide 
a non-Markovian reservoir that generates additional, de- 
layed nonlinearity, as well as spontaneous and thermal 
noise. The coupling is based on the experimentally de- 
termined Raman gain a^{uj) |48j . 

The simulation of pulse propagation entails the solu- 
tion of time-domain dynamical evolution in a quantum 
field theory with large numbers of interacting particles. 
We achieve this here primarily with a truncated Wigner 
technique 49] , which provides an accurate simulation of 
the quantum dynamics for short propagation times and 
large photon number. The quantum effects enter via 
initial vacuum noise, which makes the technique ideally 
suited to squeezing calculations. We compare simulation 
and experiments to find excellent agreement over a wide 



range of pulse energies and fibre lengths. From the simu- 
lations, we can identify the particular noise sources that 
are the limiting factors at high and low input energy. 

We begin in Sec. |TT] with an introduction to polarisa- 
tion squeezing by means of the Kerr effect, from a single- 
mode picture, before presenting the detailed model of 
pulse propagation in fibres in Sec. IIIII Sections IIVI and 
|V] describe the numerical simulation methods used, while 
the experimental set-up is described in Sec. |Vll Sec. IVIII 
discusses the results of both the experiment and simu- 
lations. The appendices contain further details of the 
theoretical description and numerical simulation. 



II. SQUEEZING 

A. Kerr squeezing 

The generation of squeezed optical beams requires a 
nonlinear interaction to transform the statistics of the in- 
put, which is typically a coherent state. The first observa- 
tion of quantum noise squeezing used four wave mixing [3] 
arising from the third order electric susceptibility x ■ 
Although material dispersion can place limits on the in- 
teraction length, this limitation can be circumvented by 
use of degenerate frequencies |50|. as in the optical Kerr 
effect. Here the interaction has the effect of introducing 
an intensity dependence to the medium's refractive index, 
Eq. ([1]), which in turn induces an intensity-dependent 
phase shift in incident pulses. This effect dominates the 
nonlinearity in fibers made of fused silica, a material with 
an inversion symmetric molecular ordering. In a pure 
Kerr material the refractive index is an instantaneous 
function of the optical intensity and the refractive index 
n is then given to second order by [5l|: 
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(1) 



where the optical intensity is given by / = ^noeoc\E\'^ 

(3) 

and Xxxxx is the third order susceptibility coefficient for 
the degenerate mode x. The instantaneity of fused sil- 
ica's nonlinearity is true only to a first approximation. 
In reality, it is only the electronic contribution, which 
typically comprises 85% of the total nonlinearity 52^, 
that is instantaneous on the scale of the 130 fs pulses 
used here. The time dependence of the remainder can- 
not be neglected and arises primarily from Raman scat- 
tering 53\. Nonetheless, the simplification of an instan- 
taneous response can be useful in gaining physical insight 
into the Kerr squeezing mechanism. 

Figure [ija) illustrates the effect of an instantaneous 
nonlinear refraction. Sending an ensemble of identical co- 
herent states into a perfect Kerr medium causes a distor- 
tion of the initially symmetric phase-space distribution. 
One can explain this distortion by considering the input 
to consist of a superposition of photon number states, 



which the Kerr effect rotates relative to one another in 
phase space. The initially symmetric phase-space dis- 
tribution characteristic of coherent states is thereby dis- 
torted into an ellipse or 'squeezed' circle. Generally the 
squeezed state will be crescent shaped, however for the 
experimental conditions of high intensities and small non- 
linearities our states never become significantly curved. 
The resultant quantum state is quadrature squeezed, 
where the squeezed quadrature X{9sq) is rotated by 6sq 
relative to the amplitude quadrature or radial direction. 
The state's phase-space uncertainty distribution is al- 
tered such that the statistics in the amplitude quadra- 
ture remain constant in keeping with energy conserva- 
tion. Thus the squeezed or noise-reduced optical quadra- 
ture cannot be detected directly in amplitude or intensity 
measurements. A detection scheme sensitive to the angle 
of the squeezed ellipse 9sq is required. 
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FIG. 1: (a) Representation in phase space of the evolution 
of a coherent beam (bottom right) under effect of the Kerr 
nonhnearity, which generates a quadrature (or Kerr) squeezed 
state (upper left). The arrow indicates the direction of state 
evolution with propagation, (b) Polarization squeezing gener- 
ated by overlapping two orthogonally (i.e. x- and y-) polarized 
quadrature- squeezed states. 



B. Single- mode picture of polarisation squeezing 



where i,j, A: = 1,2,3 and where e is the antisymmetric 
symbol. These commutation relations lead to Heisen- 
berg inequalities and therefore to the presence of intrinsic 
quantum uncertainties in analog to those of the quadra- 
ture variables. However, the fundamental noise limit de- 
pends on the mean polarization state: 



A^S'.A^S', > e 



ijk 



(Sk) 



(4) 



where the variance of Si is given by A'^Si ~ (Sf) — {Si)'^. 
This quantum picture of the polarization state of light 
cannot be represented as a point on the Poincare sphere, 
but rather as a distribution in the space spanned by the 
Poincare parameters, analogous to the phase-space rep- 
resentation of quantum optical states. This is visualized 
in Fig. [21 which shows the variances, i.e. full-width at 
half-maximum of the marginal distributions, of a coher- 
ent and a polarization squeezed state. 
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FIG. 2: Representation of the variances of a polarization 
squeezed (upper left) and a coherent state (lower right) on 
the Poincare sphere. 



The characterization of quantum polarization states 
relies on the measurement of the quantum Stokes opera- 
tors (see Rcf. 35] and references therein). These Hermi- 
tian operators are defined analogously to their classical 
counterparts as |54| : 



^n 



= alax + a^yCiy, 



Si = a^ax — ala 



y"-y^ 



S2 = alciy -\- alax , S3 ^ i{alax - alliy) , (2) 

where Qx and dy are two orthogonally polarized modes 
(with temporal, position and mode dependence implicit). 
These operators obey the SU(2) Lie algebra and thus, 
within a factor of 2, coincide with the angular momen- 
tum operators. The commutators of these operators, fol- 
lowing from the noncommutation of the photon opera- 
tors, are given by: 



Sq, Si — and 



Si,Sj — Zieijf^bk, 



(3) 



Despite the fact that the Stokes uncertainty relations 
are state dependent, it is always possible to find pairs 
of maximally conjugate operators. This is equivalent to 
defining a Stokes basis in which only one parameter has 
a nonzero expectation value. This is justified insomuch 
that polarization transformations are unitary. Consider 
a polarization state described by (Si) — (Sj) = and 
{Sk) = (So) "^ where «, j, k represent orthogonal Stokes 
operators. The only nontrivial Heisenberg inequality 
then reads: 



A^^.A^^, >|(5fe)P = |(5o)|' 



(5) 



which mirrors the quadrature uncertainty relation, and 
polarization squeezing can then be similarly defined: 



A^^, <|(5fc)|<A2^, 



(6) 



The definition of the conjugate operators Si,Sj is not 
unique and there exists an infinite set of conjugate op- 
erators S±{9),S±{6 + |) that are perpendicular to the 

state's classical excitation Sk, for which (S±{9)) — for 
all 9. All these operator pairs exist in the Si — Sj 'dark 
plane,' i.e. the plane of zero mean intensity. A general 
dark plane operator is described by: 



S±{9) = cos{9)S.,+sini9)Sj, 



(7) 



where 9 is an angle in this plane defined relative to Si 
Polarization squeezing is then generally given by: 



A^S^{9sg) < \{So)\ < A^S^{9sg + -), 



(8) 



where S±{9sq) is the maximally squeezed parameter and 
Sj_{9gq + f) the antisqueezed parameter. 

Consider, for example, the specific case of a +S3 or 
cr+-polarized beam as in the experiments presented here. 
Let this beam be composed of the two independent modes 
Qx , CLy with a relative ^ phase shift between their mean 
values. This is depicted in Fig. [T]Jb) and described by 
{cLy) = i{ax) = ia/V2 and a £ R. The beam is then 
circularly polarized with a^^ as the mean field and aa_ 
is the orthogonal vacuum mode: 

acr+ = ^ {cix - ia,y) with (acr+) = -a, 
v2 

aa_ = —={ax+iay) with {a„+)=0- (9) 

v2 

The Stokes operators in the plane spanned by Si — S2 
correspond to the quadrature operators of the dark dcr_ - 
polarization mode. Assuming \{Sa)\ ^ |a| and consider- 
ing only the noise terms, we find: 

68^(9) = a(^5a„_e-''' + Sai_e'''^ 
= aSX„_{9) 



= a {SXx{9) + 6Xyi9 - -)) , (10) 



where the Stokes operator definitions of Eq. ([2]) have been 
used in a linearized form. The sum signal, a measure of 
the total intensity, is given by: 



5S 



= a 



(^ScLa^+Sai^'^ :^a6X^^, (11) 



and thus exhibits no dependence on the dark mode. This 
considering of the physical interpretation of polarization 
squeezing shows that polarization squeezing is equiva- 
lent to vacuum squeezing in the orthogonal polarization 
mode: 

A^S^{9) < |(a)|2 ^ A^X,_{9) < 1. (12) 

Whilst a particular case is considered here, a straight- 
forward generalization to all other polarization bases is 



readily made as polarization transformations are unitary 
rotations in SU(2) space. 

In dark-plane Stokes measurements, the beam's inten- 
sity is divided equally between two photodetectors. Such 
measurements are then identical to balanced homodync 
detection: the classical excitation is a local oscillator for 
the orthogonally polarized dark mode. The phase be- 
tween these modes is varied by rotating the Stokes mea- 
surement through the dark plane, allowing full charac- 
terization of the noise properties of the dark, y-polarized 
mode. This is a unique feature of polarization measure- 
ments and has been used to great advantage in many 
experiments, for example [33, [33, [33, [5^, [Sfl [53, [H, \Em ■ 
This has also allowed the first characterization of a bright 
Kerr squeezed state as well as the reconstruction of the 
polarization variable Q function using polarization mea- 
surements [60( 1 ■ 

To show how an S3 polarised state is squeezed by the 
Kerr effect, we consider the essential Kerr Hamiltonian: 



t;; \'2 j.{nU ^2 



H = (alax) + {alay) 



(13) 



which in terms of the Stokes operators, can be expressed 
as 



H 



\{sf. + si} 



(14) 



The first term is a constant of the motion, since So gives 
the total number of photons, and has no effect on the 
dynamics. The second term is a nonlinear precession 
around the Si axis: the rate of precession is proportional 
to Si , which is a manifestation of the intensity-dependent 
refractive index of the Kerr effect. The nonlinear preces- 
sion will distort an initially symmetric distribution cen- 
tred in the iS'i-iS'2 plane (eg the 6*3 circularly polarised 
state located at a pole of the sphere) into an ellipse. 
As for ordinary quadrature squeezing, the nonlinear pre- 
cession preserves the width in the ^i direction, and so 
the squeezing is not directly observable by a number- 
difference observation. 

The advantage of the squeezed 5*3 state, as opposed to 
squeezing of a linearly polarised S2 state, is that a sim- 
ple rotation around the S3 axis allows the squeezed (or 
antisqueezed) axis of the ellipse to be aligned to the Si 
axis and thus to be detected with a number-difference 
measurement. Such a rotation is easily implemented ex- 
perimentally with a polarisation rotator. 



III. PULSE PROPAGATION 

A. Multimode description 

We have so far described the polarisation squeezing as 
a single-mode Kerr effect. However, this is accurate only 
for CW radiation, corresponding to a single momentum 
component. Ultrashort pulses, on the other hand, corre- 
spond to a superposition of many plane waves and thus 



require a multimode description. Such a description is 
crucial for an accurate treatment of dispersive and Ra- 
man effects. For a continuum of momentum modes, we 
can express the superposition as 



*,(i,z) 



1 



dfca<,(t,fc)e*('=-'=«)^+*'^°*, (15) 



where instead of annihilation and creation operators for 
each polarisation mode, we now have field operators 
\I'J.(i, z), \I'o-(i, z) for the envelopes of each of the po- 
larisation modes a = {x,y). The commutation relations 
of the fields are 



^,{t,z),^l{t,z')] ^S{z~z')S, 



(16) 



and with this normalisation, 
ber of cr-photons in the fibre 

rL 



the total nuni- 
is thus Nait) = 
J^''dz'^l{t,zy^,{t,z). 

The general quantum model for a fibre with a single 
transverse mode is derived in 6l|. The relevant aspects 
for the current system include the dispersive pulse prop- 
agation, the electronic polarisation response that gives 
the instantaneous x i ^-nd the nonresonant coupling to 
phonons in the silica. 



B. Electromagnetic Hamiltonian 

In terms of the field operators for the slowly varying en- 
velope defined above, the normally ordered Hamiltonian 
for an electromagnetic pulse in a polarisation-preserving 
fibre under the rotating-wave approximation is: 



Hem - hJ2 dz dz'uj^iz ~ z')%l{t,z)%„{t,z') 

-hxEj2jdz^i\t,z)^lit,z), (17) 



where to^lz) is the Fourier transform of the dispersion 
relation: 



OJ, 



„{z) = l- /dfc^,(fc)e^('=-'=«)^ 
Zn I 



(18) 



and xe is the strength of the third-order polarisation re- 
sponse. The birefringence of the polarisation response 
means that there are differences between the dispersion 
relations cux and ujy. The x term is assumed to be 
independent of polarisation, and cross-Kerr effects are 
neglected, as the different group velocities of the pulses 
mean that the length of time that the pulses overlap in 
the fibre is negligible. The fibre is assumed to be ho- 
mogeneous, with both uja{k) and xe independent of the 
distance z down the fibre. 

To simplify the description of the dispersive part, we 
Taylor expand uj(k) around fc = fco up to second or- 
der, which introduces the group velocity v = duj / dk\k=ko 



and dispersion parameter w" = d'^u}/dk'^\k=ko- Sub- 
tracting off the free evolution at the carrier frequency 
LjQ = ujx{ko), one obtains the simplified Hamiltonian: 

^EM = hJ2jdz^^iV^^,-^V^,) 

+ ^v$tv$. - hxE^i^^l'^ . (19) 

Here we have not included the difference in phase veloc- 
ity between the two polarisations, which just leads to a 
constant relative phase shift. 

For the methods that we use in this paper, it is conve- 
nient to treat the quantum dynamics in the Heisenberg 
picture, with time-evolving field operators. The equation 
of motion of the field annihilation operator that arises 
from the electromagnetic Hamiltonian is 



dt 



-*n 



*.,i?E 



EM 



= \-VaV+^-V^+lXE^a^i}^a. (20) 



C. Raman Hamiltonian 

As well as the interaction with electrons that produces 
the polarisation response, the radiation field also inter- 
acts with phonons in the silica. The photons can ex- 
cite both localised oscillations of the atoms around their 
equilibrium positions (Raman effect) as well as guided 
acoustic waves (GAWBs) along the waveguide. The lat- 
ter can be treated as a low-frequency component of the 
Raman spectrum, and produces random fluctuations in 
the refractive index. However the effect of this is largely 
removed in this experiment through common-mode re- 
jection, and any residual phase-noise can be accounted 
for by simple scahng laws (see section FVHCp . 

The Raman interactions produce both excess phase 
noise and an additional nonlinearity. The atomic oscilla- 
tion is modelled as a set of harmonic oscillators at each 
point in the fibre, and is coupled to the radiation field by 
a simple dispersive interaction: 

Hn = hY^Rk J dz%U^.)%Az) {bak{z) +V^,{z)} 

(7,k 

+hJ2^kbU^)b.k{z), (21) 

a.k 

where the phonon operators have the commutation rela- 
tions 

b^k{z,t),bl^,{z',t)\ = <5(z-z')4,fc'<5.,.' (22) 

The spectral profile of this interaction R{uj) is well 
known from experimental measurements [53| and is sam- 
pled here by oscillators of equal spectral spacing Auj = 



(jJk+1 — ijJk, such that hniAtj^o -Rfc/vAw — R{LUk). The 
finite spectral width of the Raman profile means that the 
Raman contribution to the nonlinearity is not instanta- 
neous on the time-scale of the optical pulse, leading to 
such effects as the soliton frequency shift |62l . |63| . 

With the Raman and electromagnetic Hamiltonians 
combined, one can derive complete Heisenberg operator 
equations of motion for the optical field operator and the 
phonon operators [6l|: 



The Wigner method, on the other hand, is an ap- 
proximate method that is valid for large photon num- 
ber n and short fibre length L. Here it is symmetrically 
ordered correlations that correspond to stochastic aver- 
ages. Because of this symmetric ordering, the quantum 
effects enter via vacuum noise in the initial conditions [7a|. 
making it a simple and efficient method for squeezing 
calculations 49] . It is also enjoying increasing utility in 
BEC simulations I79I. 






;E^4^^ + ^^] 



■" 92 



dt 



bkiz,t) 



CLU U TP " -i- " 

2 az^ 
-iLUkbk-iRk'^^i', (23) 



where we have suppressed the polarisation index, since 
the equations for each polarisation are independent. 
The initial state of the phonons is thermal, with 

(Si,(z',0)Sfc(z,0)) =nthK)4,fe"5(^-^'), (24) 

where rith(^) — 1/ [exp{hu! / kT) — 1] is the Bose-Einstein 
distribution. 



IV. SIMULATION METHODS 
A. Phase-space methods 

Phase-space methods are a means of simulating the 
dynamics of multimode many-body quantum systems. 
They are based on (quasi)probabilistic representations of 
the density matrix that are defined by means of coherent 
states. Because they are based on coherent states, they 
are ideally suited to simulating quantum optical experi- 
ments, which in so many cases begin with the coherent 
output of a laser. The two representations th at g ive rise 
to practical numerical methods are the +P [6l . l65l .l66l.l67t 
and Wigner 168* distributions. In both methods, the re- 
sultant description has the same structure as the mean- 
field, or classical, description, which is a form of nonlinear 
Schrodinger equation in the case of optical fibres. How- 
ever there are also additional quantum noise terms, which 
may appear in the initial conditions or in the dynamical 
equations. 

The +P method provides an exact probabilistic de- 
scription in which stochastic averages correspond to nor- 
mally ordered correlations. Because of this normal order- 
ing, it is suited to intensity correlation measurements. 
Quantum effects enter by stochastic terms that have 
the form of spontaneous scattering. The +P method 
has been applied to a variety of quantum-optical appli- 
cations, including superfluorescen ce 1691 |70|, parametric 
amplifiers 71] and optical fibres [14|.|47|. More recently, it 
has been applied to a variety of Bose-Einstein condensate 
(BEC) simulationsfzllzi&lziSIzS] 



B. Wigner equations 

The Wigner representation maps the operator equa- 
tions of motion onto (almost) equivalent stochastic 
phase-space equations. The mapping is not exact be- 
cause the 'nonlinear' term leads to higher-order (higher 
than second) derivatives in the equation for the Wigner 
function, which must be neglected in order that the map- 
ping to stochastic equations can be completed. These 
neglected terms are the ones, for instance, which would 
allow the Wigner function to become negative. 

The resultant equations are, up to a constant phase 
rotation, 



l*(^'^) = 



d 



-^E^^-i^'^' + ^fci-^'T^ 



dz 



2 dz 



zX^Jv!/] 



* 



dt 



bkiz,t) 



-iujkbk - iRk I* 



2A2 



(25) 



where we have assumed that the fields will be discretized 
over a lattice with segment size Az. The initial condi- 
tions are 



6fe(z,0) - rt(z) 

vi'(z,o) = (#(z,o)) + r*(z). 



(26) 



where the stochastic terms have correlations 



(rt/(z')r^(z)) = \n,^{u^k) + ^\6k,k'5{z-z') 



(r^(z')r*(z)) = 



5{z-z') 



(27) 



C. +P equations 

Phase-space equations that correspond exactly to the 
operator equations can be defined over a doubled phase 
space using the +P representation. Quantum effects en- 
ter here through multiplicative noise terms in the equa- 
tions, which generally lead to a larger sampling error than 
the Wigner method for squeezing calculations. While the 



Wigner method was used for nearly all of the simula- 
tions presented here, the +P method provides important 
benchmark results, and was used to check the validity of 
the Wigner calculations in key cases. 
The resultant +P equations are 



l*(^'*) = 



k 



d 



Rk{bk + b+}Au;~v— 



2 az^ 



iT' 



dt 



bk{z,t) 



-iujkhk - iRk"^^"^ 



TfcS 



* 



(28) 



with equations for "^^ and 6^ that have a conjugate form 
but with some independent noise terms. The initial con- 
ditions are 



6fe(z,0) 
^(2,0) 






(29) 



where the stochastic terms have correlations 



{Ti,^ {z')Tl{z)) = nti.{ujk)5k±'5{z-z') 
{T^{z',t')T^{z,t)) = x''5{z-z')6{t-t') 

= {T^+(z',t')T^+{z,t)), 

(r«(z',i')rf(^,i)) = Rk5{z - z')5[t - 1') 

= (r«+(z',t')rr(z,t)), (30) 

with all other correlations zero. 

In writing down explicit equations for the phonon vari- 
ables, we have followed the approach of CarterjT^- In 
this approach there is some freedom in how the Raman 
noise is distributed between the photon and phonon vari- 
ables, a fact which could be exploited to optimise the per- 
formance of the simulations. The alternative approach, 
as in J61|, analytically integrates the phonon variables 
out, to give nonlocal equations for the photon fields. 



D. Scaling 



For ti^tg ^ Zg, the rescaled Wigner equations are: 



l^(Cr) 



d 



~i ^ {/3fc exp(-iJ7r) + [31 exp(if7r)} Af] 

fc 

'2 9r2 






,/.(C,r)^-^.^(^|^|^-^^jexp(.^r), (31) 

with initial conditions 

/3fe(C,r = -oo) = rf(C) 
0(C = O,t) = .fe.I/(0,ioT))+r^(r), (32) 



where the stochastic terms have correlations 



(rr(C')rf(C)) 
(r;(r')r,(r)) 



^ \ nk 
n 



Sk,k' 



s{C-C) 



S{t 



2n 



(33) 



where Uk = nthikAfl/to) 

For numerical convenience, the field is split into two 
parts - a coherent field that obeys the classical equations 
of motion, and a difference field that contains the stochas- 
tic evolutions. The equations of motion of each part are 
given in appendix I A 21 

The rescaled +P equations follow similarly from Eqs 
(f28l) -([30 |) . and are given in appendix lA 31 Because of the 
much larger sampling error that arises in the +P cal- 
culations, we make use of the fact the Wigner method 
calculates the linearised evolution exactly, and use the 
+P method only to calculate the difference between the 
linearised and full evolution. If (j)wL is a Wigner solution 
to the linearised equations, and (/)pL and (/)p are +P so- 
lutions to the linearised and full equations, respectively, 
calculated with identical noise sources, then the final so- 
lution is (j) = (j)p — 4>pL + 4>WL- Because the difference 
between the full and linearised solutions is small, (j)p and 
(j)pL have very similar fluctuations in a given run; taking 
the difference removes most of the large +P fluctuations, 
and adds in only the small Wigner fluctuations. 



To simplify the numerical calculation, we transform 
to a propagating frame of reference with dimension- 
less variables: r = (i — z/v)/tQ, C = z/zo and Vl — 
ujto, where zq = to/k". The fields are also rescaled: 
4> — '^y^vto/n and /3fe = r^bk exp(iOr)-\/zo/io'^, where 
Tk — Rky/nzo/tov'^ is the rescaled Raman coupling, 
which is related to the Raman gain a^{n) via rk = 
^Ja^\kASl)/2TT. The quantity n — v^Iq/xzo gives the 
typical number of photons in a soliton of width tg- The 
effective nonlinearity that gives rise to solitons has both 
electronic and Raman contributions: x = Xb+X-Rj where 
the Raman contribution is estimated to be a fraction 
/ = xIxr — 0-15 of the total. 



V. OUTPUTS AND MOMENTS 

We find that good precision (a few percent of the 
squeezing in decibels) is obtained when averages are cal- 
culated using 1000 realisations of the Wigner equations. 
For further precision, 10,000 trajectories can be used, in 
which case we find that the sampling error cannot be dis- 
tinctly plotted on the graphs. With the +P method, on 
the other hand, we find that at least 10,000 trajectories 
are needed in some cases to produce useful results, even 
when the differencing method is used. 

The observable moments in the polarisation squeez- 
ing measurements are integrated intensity measurements 



and their variances, which are neither simply normally 
ordered nor symmetrically ordered. Thus the results of 
the phase-space simulations must be adjusted for reorder- 
ing, as we describe below. 

In the theoretical description of the system, there are 
two optical fields, corresponding to the two polarisation 
modes of the fibre: '^x{t,z) and '^y{t,z). To describe 
the polarisation squeezing, we define integrated Stokes 
operators, which are a generalisation of Eq. ([2]): 

So = Nx4T) + Nyy{T), Sl = Nx.4T)~Nyy{T), 

S2 = NxyiT) + Nyx{T), S3=iNyx{T)- iN,y (r034) 

where T is the propagation time down the length of fibre 
and Naai{t) = J dz^'l{t,z)^a'it,z). After the polarisa- 
tion rotator, the fields are transformed to 

$^(i,z) = cos{e/2)^^{t,z)-ism{e/2)^y{t,z) 
^'y{t,z) = sin(6'/2)$^(t,z)-f icos(6l/2)$j^(i,z),(35) 

which leaves Sq unchanged but which transforms Si to 

Sg = cos(6i)5i + sin(6i)52. (36) 

To calculate that squeezing in Sg, we need to calculate 
the mean (Sg') and mean-square (Sg). 

A. +P Moments 

For the +P method, stochastic averages of the phase- 
space variables give normally ordered moments. Thus 
the mean (^Sg) can be calculated directly, as it already 
normally ordered. The mean square, however, requires a 
reordering: 



(^.^> 



cos(6l)5i+sin(6')S'2 



(^o), (37) 



as shown in appendix IB II 

For convenience, we define corresponding stochastic 
polarisation parameters Sj , sg in terms of the normalised 
+P fields: n^a'iO = / (iT0+(r, C)(/>a'(T, C). The mea- 
sured variance can then be written: 



vaT{Sg) 



re 



1 



(^^)+p-^«)+P + =N)+p ,(38) 



where (. . .) , „ denotes a stochastic average with respect 
to an ensemble of +P trajectories. The correction term 
here corresponds to the shot-noise level of a coherent 
state (for which {sf} p = (sg) ): var(S'e)coh = (So) = 



given by: 



Thus the amount of squeezing in decibels is 



Squeezing(dB) — 10 log- 



(4) +p - ^{se)lp + (sq) 



+p 



{'^o)_ 



(39) 



B. Wigner Moments 

Stochastic averages in the Wigner method correspond 
to symmetrically ordered products, thus making a re- 
ordering necessary for both the mean and variance of 
the integrated intensity measurements. First we note the 
symmetric form of Na-a' '■ 



N„ 



^ fdz[^l{z)^,iz) + ^^,iz)^U^)} 



N„ 



1 



(40) 



where M is the number of Fourier modes used to decom- 
pose the pulse shape. Because 52 and 6*3 contain only 
cross-polarisation coherences, there is no correction from 
reordering. In Si, the corrections from horizontal and 
vertically polarised terms cancel out. Thus it is only the 
total intensity that requires a correction, and this corre- 
sponds to the vacuum-energy contribution: 



So = So + M 

syni 



(41) 



The variance of the Stokes operators contain terms 
with products of four operators, each of correspond to 
24 possible orderings. As appendix IB 21 shows, most of 
the corrections cancel out, leaving: 



Sa 



= Si 



\"^ 



(42) 



Similarly to above, we can define an analogous stochas- 
tic polarisation parameter sg in terms of the normalised 
Wigner fields: no.o-'(C) = /drt/)* (r, C)0<t'(t, C)- The mea- 
sured variance can then be written: 

var(S,) = if (^{sl)^ {sg)l, - ^M^ , (43) 

where (. . .\.^ denotes an stochastic average with respect 
to an ensemble of Wigner trajectories. The shot-noise 
reference level is given by var(S'e)coh = {So) = ^{so)^/ — 
M . Thus the amount of squeezing in decibels is 



Squeezing (dB) = 10 log ■ 



-<sl) 



w 



Ww-^M 



{so\ 



LM 



(44) 



VI. EXPERIMENT 



The laser system used in these experiments is a home- 
made solid state laser where Cr^+:YAG is the active 
medium [80] ■ This system emits pulses with temporal 
widths of To =130-150 fs at a central wavelength Aq be- 
tween 1495-1500 nm. These ultrashort pulses exhibit a 
bandwidth limited secant-hyperbolic spatial amplitude 
envelope and are thus assumed to be solitons. The laser 
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FIG. 3: Schematic of the single-pass method for the efficient 
production of polarization squeezed states. The Stokes mea- 
surement after the fiber scans the dark S1-S2 plane of the 
circularly {S3) 7^ polarized output. 



repetition rate is 163 MHz and the average output power 
lies between 60 and 90 mW corresponding to pulse ener- 
gies of 370-550 pJ. 

In the present configuration, pictured in Fig. [3l laser 
pulses are coupled into only one end of the glass fiber. 
This produces quadrature squeezing rather than am- 
plitude squeezing which is not directly detectable (see 
Fig.[T]Ja)). However, overlapping two such independently 
and simultaneously squeezed pulses after the fiber allows 
access to this quadrature squeezing by measurement of 
the Stokes parameters (Fig.[ljb)). This requires the com- 
pensation of the fiber birefringence, which we choose to 
carry out before the fiber to avoid unnecessary losses to 
the squeezed beams. The optical fiber used was the FS- 
PM-7811 fiber from 3M, chosen for its high birefringence, 
i.e. good polarization maintenance, as well as its rela- 
tively small mode field diameter, i.e. high effective non- 
linearity and thus low soliton energy. The most relevant 
fiber parameters are listed in Table [H 



Parameter 


Symbol 


Fibre I 


Fibre II 


Units 


Mode field diameter 


d 


5.42 


5.69 


/im 


Nonlinear refractive 
index (xlO"^° ) 


n2 


2.9 


2.9 


m VW 


Efltective nonlinearity 
(xlO-^) 


7 


5.3 


4.8 


l/(m-W) 


Soliton energy 


Esoi 


56 


60 


pJ 


Dispersion 


k"{^P2) 


-10.5 


-11.1 


fs /mm 


Attenuation @ 1550 nm 


a 


1.82 


2.03 


dB/km 


Beat length 


Lb 


1.67 


1.67 


mm 


Polarization crosstalk 
per 100m 


AP 


< -23 


< -23 


dB 



TABLE I: Values for the material parameters for the 3M FS- 
PM-7811 fiber. Fibres I and II refer to two different pro- 
duction runs. All values (when not otherwise stated) are for 
Ao = 1500 nm and tq = 130 fs. 

For experimental ease, the polarization of the beam af- 



ter the fiber was set to be circular, e.g. a+. The orthog- 
onal Stokes parameters in the dark S1-S2 plane, given by 
Eq. [3 are measured by rotating a half-wave plate before 
a polarising beam splitter, as in Fig. [31 Equations (fTUl ^ 
llip provide an interpretation of the classical excitation in 
dcr^ as a perfectly matched local oscillator for the orthog- 
onally polarized dark mode &„_ . The phase between Oq-^ 
and d^_ varies with the rotation of the half-wave plate 
angle, 8, to give the phase-space angle 9 — 40. This 
noise level was compared with the respective Heisenberg 
limit. The sum photon current. So, gives the amplitude 
noise of the input beam, for a Kerr-squeezed state this 
equals the shot noise. This reference level was verified 
by observation of the balanced homodyne detection of a 
coherent state as well the sum of the balanced homodyne 
detection of the x- and j/-polarized modes. 

The polarizing beam splitter outputs were detected by 
two balanced photodetectors based on pin photodiodes. 
The detectors had a DC output (<1 kHz) to monitor 
the optical power as well as an AC output (5-40 MHz). 
This frequency window was chosen to avoid low fre- 
quency technical noise and the high frequency laser rep- 
etition rate. The sum and difference of the detectors' AC 
photocurrents, representing the noise of different Stokes 
variables, were fed into a spectrum analyzer (Hewlett- 
Packard 8595E) to measure the spectral power density 
at 17.5 MHz with a resolution bandwidth of 300 kHz and 
a video bandwidth of 30 Hz. 



VII. RESULTS - EXPERIMENT AND 
SIMULATION 

A. Characterising the single-pass method 

The single-pass squeezing method allows the mea- 
surement of greater squeezing as well as the direct 
and full characterization of the bright Kerr-squeezed 
beams jSSi, ^] . Both of these traits are visible in Fig. H) 
Here the measured AC noise as a function of the rotation 
of a half- wave plate (by the angle 8) in the dark Stokes 
plane is seen. A progression between very large noise and 
squeezing is observed, as expected from the rotation of a 
fiber squeezed state. Plotted on the x-axis is the projec- 
tion angle 9, i.e. the angle by which the state has been 
rotated in phase space, inferred from the wave plate an- 
gle {6 — 48). Here pulses of 83.7 pJ were transmitted 
through 13.3 m of optical fiber and the electronic signals 
were corrected for the —86.1 ±0.1 dBm dark noise. 

For 9 = 0, an 5*1 measurement gives a noise value 
equal to the shot noise. This corresponds to the am- 
plitude quadrature of the Kerr-squeezed states emerging 
from the fiber. Rotation of the state by 9sq makes the 
state's squeezing observable by projection onto the minor 
axis of the noise ellipse. Further rotation brings a rapid 
increase in noise as the excess phase noise, composed of 
the antisqueezing and the classical thermal noise arising 
from GAWBS, becomes visible. The maximum noise is 
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FIG. 4: Noise power against phase-space rotation angle for 
the rotation of the measurement half-wave plate for a pulse 
energy of 83.7 p J using 13.3 m of fiber I. Inset: Schematic of 
the projection principle for angle 6. 



observed for 9 = Oaq + ^ . Under the assumption of sta- 
tistically identical but uncorrelated Kerr-squeezed states, 
this measurement is equivalent to the characterization of 
the individual squeezed states using standard local oscil- 
lator and homodyne detection methods. However here no 
stabilization is needed after production of the polariza- 
tion squeezed state. This is advantageous for experiments 
with long acquisition times, i.e. state tomography, and 
has indeed allowed the reconstruction of the Wigner func- 
tion of the dark Stokes plane or Kerr-squeezed states [60| ■ 




0.2 0.4 0.6 0.6 

Optical transmision 

FIG. 5: Linear noise reduction against optical transmission 
for the polarization squeezing generated by pulses of an energy 
of 81 pJ in 3.9 m of fiber I. 

It is crucial to ensure that the measured squeezing did 
not arise from detector saturation or any other spurious 
effect. This was accomplished observing the noise of a 
variably attenuated squeezed beam, where a plot of the 
linear relative noise against the transmission should be 
linear for true squeezing. A representative plot for a 81 pJ 
pulse in a 3.9 m fiber exhibiting — 3.9±0.3 dB of squeezing 
is shown in Fig. [51 the linear result is indicative of genuine 
squeezing. 

The single-pass polarization squeezer exhibits a good 
temporal stability, highlighted by the results in Fig. [S] 
Here the sum (shot noise) and difference (polarization 




40 60 

Time (minutes) 

FIG. 6: Plot showing a stable squeezing o/— 4.0 ±0.3 dB over 
100 minutes. A 30 m optical fiber with a pulse energy of 80 pJ 
was used. 



squeezing) channels have been plotted. An average of - 
4.0 dB of squeezing corrected for —8 5. 8 ±0.1 dBm of dark 
noise was measured over 100 minutes. The squeezer used 
30 m of optical fiber into which two orthogonally polar- 
ized pulses of 40 pJ each were coupled. The most sensi- 
tive factor in this setup is the locking of the birefringence 
compensator. Further important parameters are the cou- 
pling of light into the fiber and the laser power stability. 
All of these parameters are easily held stable by exploit- 
ing commercially available components. 



B. Squeezing Results 

The squeezing angle, 6sq and the squeezed and anti- 
squeezed quadratures were experimentally investigated 
as a function of pulse energy from 3.5 pJ to 178.8 pJ, as 
plotted in (a), (b) and (c), respectively, of Fig. [7] (trian- 
gles). The X-axis shows the total pulse energy, i.e. the 
sum of the two orthogonally polarized pulses comprising 
the polarization squeezed pulse. We observe maximum 
squeezing — 6.8 ± 0.3 dB at an energy of 98.6 pJ. The 
corresponding antisqueezing of this state is 29.6 ± 0.3 dB 
and the squeezing angle is 1.71°. As the optical energy 
goes beyond 98.6 pJ, the squeezing is reduced, eventually 
reaching the shot noise limit (SNL), and the increment 
of antisqueezing slows down to a plateau area. 

The loss of the set-up was found to be 13%: 5% from 
the fibre end, 4.6% from optical clments and from the 
fibre attenuation (2.03 dB/km at 1550 nm), 2% from 
incomplete intereference between the two polarisation 
modes ( 99% visibility was measured), and 2% from 
the photodiodes. Thus we infer a maximum polarisa- 
tion squeezing of — 10.4 ± 0.8dB. The irnprovement in 
squeezing over previous implementations [33], |46[ of the 
single-pass scheme is largely due to the low loss achieved 
here. 

The theoretical simulations for the squeezing, anti- 
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squeezing and squeezing angle at different input ener- 
gies are also given in Fig. [7] by solid and dashed lines. 
As described in further detail below, the effect of ex- 
cess phase noise, such as GAWBS 46', is included by a 
single-parameter fit between the simulation and experi- 
mental data for squeezing angle (shown by solid line in 
Fig. [7]^ a)). The theoretical results for squeezing and an- 
tisqueezing then show very good agreement with the ex- 
perimental data, and are consistent with the measured 
linear losses of 13%. From the simulations, the effect of 
the GAWBS is seen to be a reduction in squeezing for 
lower energy pulses; above about lOOpJ, it has virtually 
no effect on the squeezing. Although some deviations ap- 
pear at very high input energy, the simulations also show 
the same deterioration of squeezing for higher pulse en- 
ergies as is seen in the experimental results; this effect 
does not occur in the simulations if Raman terms are 
neglected, as we discuss below. 



C. Phase-noise and GAWBS 

Excess phase noise, caused for example by depolarising 
GAWBS in the fibre, is determined for each fibre length 
by a single-parameter fit of the experimental and simu- 
lation squeezing angles. We model this by independent 
random fluctuations in the refractive index at each point 
along the fibre length. The cumulative effect on each 
pulse at a given propagation length is a random phase 
shift whose variance is proportional to the time-width of 
the pulse: 



r,C)=0o(r,C)e"', 



(45) 



where (rf^ ex io. 

Such phase fluctuations do not affect the number dif- 
ference measurement 5*1 , but they do lead to fluctuations 
in 5*2 



S-2-iS: 



2Tjn \Mr,0rdT^7jE (46) 



where rj is now taken to describe the relative (depolar- 
ising) relative phase shifts. Thus the variance relative 
to shot noise of 52 caused by phase fluctuations scales 
linearly with the energy of the pulse: 



var(52 

7¥" 



{tj^)E^ 



E 



cE 



(47) 



where the constant of proportionality c is to be deter- 
mined by the fit. Here we have assumed that the pulse 
width is a constant, independent of input energy. This 
assumption is not entirely accurate, because unless the 
energy is the soliton energy for that pulse width, the 
pulse will reshape to form a soliton, thereby altering the 
pulse width. However, for short fibre lengths, this effect 
should be small, and so we neglect it in our calculations. 
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FIG. 7: Measurement results and theoretical simulations for 
13.2 m 3M FS-PM-7811 fiber (run II) as a function of pulse 
energy: (a) the squeezing angle, (h) the squeezing and (c) an- 
tisqueezing noise. Solid and dashed lines show the simulation 
results with and without additional phase noise, with linear 
losses are taken to be 13%. The shading indicates simulation 
uncertainly. The simulation result without third-order disper- 
sion is given by the dots in (b). Diamonds represent the ex- 
perimental results, with experimental uncertainty indicated by 
the error bars in the squeezing. Both the simulation and the 
experimental errors were too small to be plotted distinctly for 
the squeezing angle and antisqueezing. The measured noises 
are corrected for —85.1 ±0.1 dBm electronic noise. 



The effect of the phase noise will be to stretch the 
squeezing ellipse in the S2 direction, according to the 
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formula: 

var(S0) ^ 

— a cos 



(^o) 



Qk) + hw^{Q - Ok) + cE sm^{0), 

(48) 

where 9k is the predicted angle from the Kerr-only 
squeezing, a is the relative Kerr squeezing and b is the 
relative Kerr anti-squeezing. These parameters are calcu- 
lated by the simulation, and are a function of the input 
energy E. The value of c is determined by fitting the 
predicted angle of maximum squeezing as a function of 
E against the observed values. The predicted angle is 
obtained from the extrema of the expression in Eq. (pS)) 
and the fit is performed with a nonlinear least squares 
method. Once the value of c is determined, new values 
of squeezing and antisqueezing are calculated from eq. 



As figure [8] illustrates, the excess phase noise has a 
substantial effect, on both the squeezing angle and the 
amount of squeezing, only at low levels of squeezing. For 
highly squeezed light, the Kerr-squeezed ellipse is more 
closely aligned to the phase quadrature, and thus the 
phase-noise merely has the effect of increasing the an- 
tisqueezing. This view is confirmed by the results in 
Fig. ([7]), where the difference between the curves with 
and without the phase-noise-fitting is discernible only at 
lower input energies. 





FIG. 8: Illustration of the effect of excess phase noise on dif- 
ferent squeezing ellipses. Dashed line gives the Kerr-squeezed 
ellispse, and solid line gives the ellipse with added phase noise. 
The effect on the squeezing and the angle is less for the ellipse 
with larger Kerr squeezing (lower ellipse). 



D. Third-order Dispersion 

The comparison between theory and experiment con- 
firms the deterioration of squeezing at large pulse energy 
caused by Raman effects in the fibre. However there is 
still some residual discrepancy between theory and ex- 
periments, which could be caused by various higher-order 
effects not included in the theoretical model. We here ex- 
plore the effect of third-order dispersion, and find that it 
accounts for some of the unexplained difference at high 
energies. 



Third-order dispersion Slj arises from the rate of 
change of curvature of the dispersion. It becomes more 
important for shorter pulses or when operating near the 
zero-dispersion wavelength [82]. In the propagation equa- 
tion, it appears as an extra term in the scaled equations: 



9C 6 9r3' 



(49) 



where B^ = k'" zq/Iq is a dimensionless third-order dis- 
persion parameter. For the fibre used in the exper- 
iment, the third-order dispersion at A = 1499nm is 
k'" = 8.38 X lO-^isVm, giving B3 = 0.097. The effect of 
third-order dispersion on the pulse spectrum for various 
energies is shown in Fig. [9l where significant differences 
appear only above the soliton energy. 

Third-order dispersion does not have an observable ef- 
fect on the squeezing angles or the amount of antisqueez- 
ing, but its effect can be seen on the squeezing, as shown 
in Fig. EKb) by the difference between the solid and dot- 
dashed lines. Below the soliton energy, the third-order 
dispersion has no observable effect. It diminishes the 
amount of squeezing at around the soliton energy, and 
at higher energies it changes the rate at which squeezing 
deteriorates as a function of energy. Far above the soliton 
energy, there remains some difference between simulation 
and experiment, which indicates that other higher-order 
processes may be playing a role at these energies. Be- 
cause, in any case, the effect of third-order dispersion is 
rather small, we do not include it in the other simulation 
results shown in this paper 



E. Raman noise effects 

The Raman effect has a significant effect on the pulse 
shape and spectrum for the more intense pulses at these 
subpicosecond pulse widths. For a soliton pulse, the ef- 
fect of the Raman interactions is to induce a frequency 
shift in the soliton and hence a delay in its arrival 
time [62, [631 . For pulses above the soliton energy, the Ra- 
man interaction affects the way the pulse reshapes. With 
a purely electronic (instantaneous) nonlinearity, the pulse 
reshapes into a narrower soliton, at the same time shed- 
ding radiation that forms a low pedestal underneath the 
soliton. In the frequency domain, this results in a modu- 
lation of the pulse spectrum. As Fig. [TUl shows, with the 
Raman terms included, the reformed soliton separates 
from the pedestal, which distorts the spectrum into an 
asymmetric shape. 

For a pure Kerr effect, the squeezing is proportional to 
the intensity of the light, which in our case corresponds 
to the input energy of the pulse. However the exper- 
imental and simulation results clearly show that while 
the squeezing increases with input energy over a range 
of energies, there is a point beyond which the squeezing 
deteriorates. This deterioration is largely due to Raman 
effects, as Fig. [TT] reveals, which compares the simula- 
tions with and without Raman effects. In the latter case 
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FIG. 9: Simulation pulse-spectrum at pulse energies (a) 
l.^Es, (b) Es, (c) O.5E3 at rescaled distances 0/ ^ = (dot- 
dashed), and ^ = 25 with (solid) and without (dashed) third- 
order dispersion 



the nonlinearity is taken to be of the same magnitude as 
the former but is instantaneous. Without Raman effects, 
the squeezing does not suffer the same catastrophic re- 
duction at high energies, but it does appear to saturate 
at around the soHton energy (2x54pJ), demonstrating 
that pulse-reshaping effects are also in play. 

For L = 13.35m, the optimum energy is around 80% 
of the soliton energy. 
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FIG. 10: Simulation pulse shape (a) and spectrum (b) at 
pulse energies 1.5Es and normalised propagation length C, = 
25, with (solid) and without (dashed) Raman effects. 



F. Comparison with exact -{-P results 

Nearly all of the simulation results presented in this 
paper were calculated with the truncated Wigner phase- 
space method, because results can be obtained quickly 
and with low sampling error. However, the Wigner tech- 
nique only provides an approximation to the true quan- 
tum dynamics. While the approximation is usually a 
good one for intense optical pulses, some deviations from 
the exact result could in principle occur for long simu- 
lation times, or when highly-squeezed states are being 
produced. To test the Wigner method, we compared 
selected points with -\-P calculations, and found agree- 
ment within the statistical uncertainty. One example is 
shown in Fig. [TTl where the -\-P results are shown as the 
squares. As the error-bar indicates, the sampling error 
for the -\-P is much larger than that of the Wigner for 
the more intense pulses, even though 10 times as many 
trajectories were used for the +P calculation. Even for 
the same number of trajectories, -\-P calculation is more 
computationally exacting. This combination of greater 
computational cost per trajectory and the larger num- 
ber of trajectories required for a meaningful -\-P result is 
why the Wigner technique has been our method of choice 
for squeezing calculations. The -\-P method comes into 
its own when more exotic quantum states or fewer pho- 
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FIG. 11: Squeezing degradation at high intensity: (a) squeez- 
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bre I. Solid and dashed hnes shows the simulation resuhs with 
and without Raman effects (i.e. a purely electronic nonlinear- 
ity), respectively. The data points marked by squares are the 
results of exact -|-P calculations with error-bars indicating es- 
timated sampling error. The points at 109 pJ and 134pJ were 
calculated with 10,000 trajectories; the other 4 -\-P points 
were calculated with 1000 trajectories. Note that these re- 
sults were obtain in an experimental set-up with higher losses 
than that of Fig.[71 giving a reduced magnitude of raw squeez- 
ing. The simulations here assume 19.9% loss 



tons are involved, i.e. when the Wigner technique is not 
expected to be reliable. It is also possible that appro- 
priate difFusion[83] or drift 84] gauges may improve the 
performance of the +P calculations. 



G. Comparison for different fibre lengths 

The squeezed and antisqueezed quadratures as well as 
the squeezing angle Qsq of such polarization states were 
investigated as a function of pulse energy for different 
lengths of 3M FS-PM-7811 fiber, as shown in Figs.[Tl[ll 
The figures are organized into pairs of lengths: Fig. [T^ 
shows 3.9 and 13.3 m. Fig. [13] shows 20 and 30 m and 
Fig. [HI shows 50 and 166 m. For each length the squeez- 
ing angle (±0.3°), squeezing (±0.3 dB) and antisqueezing 
(±0.3 dB) form a column. Due to the technical limita- 



tions of the photodetectors it was not possible to measure 
above 125 pJ or 20 mW in this particular experimental 
run. The losses in this particular set-up were also larger 
than in that which gave the results in Fig. [7] 

Even though the simulations and experiment agree 
very well for the angle and the antisqueezing, some small 
discrepancies appear in the squeezing at longer fibre 
lengths. This could be caused by variation of the ma- 
terial parameters along the fibre length, or inaccuracies 
in the Raman model, which would become more promi- 
nent for longer fibres. 

Ideal Kerr squeezing should increase with propagation 
distance. However the experimental data and simulations 
show that, above 13.4m, the squeezing at a given input 
energy is largely insensitive to the length of fibre. The 
exception here is that the deterioration of squeezing due 
Raman effects starts to occur at slightly lower energies. 
Thus, the maximum squeezing available at a given fibre 
length actually decreases for longer lengths. Meanwhile 
the antisqueezing increases with propagation distance, as 
expected. 
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FIG. 13: Experiments (corrected for dark noise) and simula- 
tions (with and without fitted phase noise) of the polarization 
squeezing, antisqueezing and squeezing angle for 20 (a, c, e) 
and 30 m (b, d, f) of fiber II. The lighter data points in (c) 
are from a corrected experimental run. 



FIG. 14: Experiments (corrected for dark noise) and simula- 
tions (with and without fitted phase noise) of the polarization 
squeezing, antisqueezing and squeezing angle for 50 (a, c, e) 
and 166 m (b, d, f) of fiber II. The amount of dispersion over 
166 m makes the simulations impractical for this case. 



H. Optimal squeezing as a function of 
power /length 



Some insight can be gained from further simulations of 
squeezing as a function of fibre length, for various input 
energies, as shown in Fig. [151 This figure reveals that 
for a given input energy there is an optimum length for 
the best squeezing, reflecting the length-dependence of 
the Raman-induced deterioration revealed in the previ- 
ous plots. The best squeezing overall is obtained for a 
pulse at the soliton energy (54pJ in each pulse), which 
indicates that the reduced optimal squeezing at other 
energies is due to pulse-reshaping effects. Thus for the 
to =130fs used here, the optimum fibre length would be 
L ~ 7m, (although the improvement over 13m would 
only be a fraction of a dB). Alternatively, for a fixed fi- 
bre length, one could optimise the maximum squeezing 
by changing the pulse-width to yield a soliton at that 
point. 

Furthermore, as Fig. [12] plots the simulation results 
without adjustment for linear loss, it shows that inferred 
squeezing of over -12dB is possible. 
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FIG. 15: Simulated squeezing as a function of fibre length, 
for various total input energies: E =74.4pJ (triangles), 
E =83.7pJ (diamonds), E =93pJ (squares), E =108pJ (cir- 
cles) and E =119pJ (crosses). Linear loss and phase- noise 
have not been included. (Fibre I) 



VIII. CONCLUSION & OUTLOOK 

An excellent —6.8 ± 0.3 dB of polarization squeezing, 
the greatest measured in fibres to date, has been demon- 
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strated with the novel single-pass setup [34|. From this 
value it is possible to inferr that — 10.4±0.8 dB of squeez- 
ing was generated in the fiber. To further improve the 
measured noise reduction, losses after the fiber must be 
minimized by, for example, employing more efficient pho- 
todiodes in a minimal detection setup using highest qual- 
ity optics. We speculate that net losses of as little as 5% 
should be possible, thereby allowing the measurement of 
squeezing in excess of -8 dB. 

By exclusion of the Raman and/or the GAWBS ef- 
fects in the simulations, it is clear that that the former 
is a limiting factor for high pulse energies, whereas the 
latter is detrimental at low energies. Investigation of a 
range of fibre lengths revealed that greater squeezing is 
not achieved going beyond 13.2m. Indeed, simulations 
indicate that slightly greater squeezing may be achiev- 
able at a lower fibre length of around 7m. 

Further improvement may be possible through the use 
of photonic crystal fibers (PCF), which are novel fibers 
manufactured with specially designed_ light-guiding air- 
silica structures along their length [85|- These have al- 
ready been used in several squeezing experiments [86|, |87|, 
ISSl. 189(1 ■ and with fewer low- frequency acoustic vibrations, 
are also expected to improve squeezing results by mini- 
mizing destructive GAWBS noise [90] ■ Such an advance 
would bring fiber-produced squeezed states closer to min- 
imum uncertainty states, a desirable feature for quantum 
information applications. 
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where the loss noise has the correlations 



[Tl (t, c)rL (t , c )) = ^ . ( A2) 



The +P equations are similar, except that there is no 
stochastic terms from the loss. 

Naturally, because this loss is not a physical effect, the 
time window should be wide enough so that the loss does 
not affect the pulse itself. 



2. Split Wigner equations 



To increase the output precision in the numerical 
Wigner calculation, we split the fields into means plus 
deviations: (t> = (t> + 5(f), fij, = [3f. + 6 (3k and evolve the 
two parts separately. The simulated equations are thus: 






dc 



r'^lM'^J^kT 



= -iri i (j}A(t)* + A(t)(j) + A(t)A(t)* 
i d'^ - - 



2nzoACj ' 



A(j) 



2 9x2 



-iA<j)*(j)'^ - i {IA(J) + AI(j) + AIAc, 



A(t> + iU* + A(t)*\ (2 * A(/)0 -f Act?) 



(A3) 



where 



APPENDIX A: NUMERICAL 
IMPLEMENTATION 

1. Absorbing potentials 

The split-step Fourier method used to integrate the 
equations gives periodic boundary conditions. Thus any 
radiation shed by the pulse that reaches the edge of 
the time window during the simulations will re-enter the 
other side, eventually interfering with the original pulse 
and thereby giving possibly spurious results. To prevent 
this, we include an inhomogenous loss that absorbs any 
shed radiation before it reaches the boundary. We choose 
a (negative) gain profile of ^(t) = sin^"(7rr/2To), where 
2to is the width of the simulation window. The contri- 
bution to the Wigner equations is then: 



;^'A(r,C) 



^~9{r)cl>{T,0+TL{rX), 



(Al) 



/(C,r) - ^25R{/?,e-*^^''^}Ar! 
A/(C,t) = ^23fi{A/3fce-*^^''^}Arj (A4) 



and with initial conditions 



^(C = 0,r) . 


= v]Vsech(T 


A0(C = O,r) . 


-- r0(r) 


^,(C,r = -oo) . 


= 


i/3fc(C,'r = -oo) = 


-' rfe(C). 



loss 



(A5) 



The soliton number is defined as A^ = E/ Eg , where Es = 
2hujn is the energy of a fundamental sech soliton of width 

to. 
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3. Rescaled +P equations 

The rescaled +P equations, corresponding to the 
Wigner equations of Eq. ([55)1 are 



d 

9C' 



^E{/^^^"'" + /3,V-}^^ + i|:2 



+z(l - /)0+(/) + VtT^iC, t) + iT^iC, t) 



Ot 



Pk = r^(-^</)+0 + ^«(C,r))e^"^ 



(A6) 



with equations of conjugate form for 0+ and P'^ . The 
initial conditions are 



-r/3 



r^(C) 



n 



Pk{C,T 

0(C = O,t) 
where the stochastic terms have correlations 



(A7) 



(rf, (C')rf(C)) 
(r^(C',T')r^(C,r)) 



rlnkSk,k' 

nAfl 
1-/ 



HC-C) 



6{C - C)6{t - t') 

(r^+(C',r')r^+(C,r)), 



(r«(C',r')r«(C,r)) = -5(C-C'M(r-r') 

= (r«+(C',r')rf+(C,r)).(A8) 

Preliminary investigation of other (physically equivalent) 
ways to numerically implement the Raman noise did not 
find any improvement over the simple choice given here. 



APPENDIX B: OUTPUT MOMENTS 

As discussed in Sec. |Vl the +P and Wigner simulation 
methods give correlations that are normally and sym- 
metrically ordered, respectively. To compare with ex- 
perimental measurements of the Stokes' parameter vari- 
ances, some re-ordering is necessary. The transformation 
^ that relates Si to Sg preserves the commutation rela- 
tions. Thus to reorder the variance of a Stokes parameter 



at a general angle in the dark plane, Sg, we only need to 
consider the corrections that arise from reordering Sf. 

1. Normal Ordering 

First we expand the mean-square of Si as 

■.Sf: = : N^^N.^^ : -2N^^Nyy+ : NyyNyy : .(Bl) 
Thus we need only consider the reordering of terms of 
the form: 



= jdzj dz'^Ut, ^)*U<, Z')^a{t, Z)^,{t, Z') 

+ I dzdl{t,z)^„{t,z) 
= : N,,N,, : +N,,, (B2) 

which gives 

Si = -.Si -.+30. (B3) 

2. Symmetric ordering 

To symmetrically order the relevant products of 4 field 
operators, we must start from the sum of all 4! = 24 
possible orderings. The essential result that we require 



sym 



(B4) 



where M is the number of modes. The mean square of 
Si is then: 



Sf 



- [N^XN^X- 2N^^Nyy + Nyy Nyy 



Sf + -M. 



(B5) 
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